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Abstract 



We report on a calculation of Bk with domain wall fermion action in 
quenched QCD. Simulations are made with a renormalization group improved 
gauge action at /3 = 2.6 and 2.9 corresponding to a^^ k, 2GeV and 3GeV. 
Effects due to finite fifth dimensional size and finite spatial size A'^,- are 
examined in detail. Matching to the continuum operator is made pertur- 
batively at one loop order. We obtain BxilJi = 2GeV) = 0.5746(61)(191), 
where the first error is statistical and the second error represents an estimate 
of scaling violation and O(a^) errors in the renormalization factor added in 
quadrature, as an estimate of the continuum value in the MS scheme with 
naive dimensional regularization. This value is consistent, albeit somewhat 
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small, with BxifJ- = 2GeV) = 0.628(42) obtained by the JLQCD Collabora- 
tion using the Kogut-Susskind quark action. Results for light quark masses 
are also reported. 
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I. INTRODUCTION 



The kaon B parameter Bx is an important quantity to pin down the Cabibbo-Kobayashi- 
Maskawa matrix from experiment, thereby advancing our understanding of CP violation in 
the Standard Model [||. 

A crucial ingredient in a precision calculation of Bk is chiral symmetry. Without this 
symmetry the relevant AS = 2 four-quark operator mixes with other operators of different 
chiralities. It is a non-trivial task to accurately determine the mixing coefficients. 

This problem has caused significant difficulties with calculations using the Wilson-type 
fermion action, which has explicit chiral symmetry breaking. While several non-perturbative 
methods have been developed to determine the mixing coefficients the numerical errors 
in the values of Bk obtained with these methods are still quite large [|I| . 

The situation is better with the Kogut-Susskind fermion action for which U{1) subgroup 
of chiral symmetry, valid at finite lattice spacings, ensures the correct chiral behavior of the 
matrix element [Q. Exploiting this feature a systematic and extensive set of simulations 
have been carried out 0. Taking into account 0{a^) scaling violation and 0{ay^) errors 
that arise with the use of one- loop perturbative renormalization factors, BxifJ' = 2GeV) = 
0.628(42) has been obtained in the continuum limit in the MS scheme with naive dimensional 
regularization (NDR). 

Recent development of the domain wall and overlap ||lO|,|lll fermion formalisms has 
opened a prospect toward an even better calculation. Even at finite lattice spacings, these 
formulations maintain both fiavor and chiral symmetries, either of which is broken in the 
Wilson-type fermion action or the KS quark action. Hence one expects that systematic as 
well as statistical uncertainties are better controlled in these formulations than others. A 
pioneering calculation of the Bk parameter in this direction was made in Ref. |jl2| using the 
domain wall fermion formalism of QCD (DWQCD). In this article we present results of our 
study toward a precision determination of Bk with DWQCD. 

Our investigation is carried out in the quenched approximation using Shamir's formu- 
lation of domain wall fermion for quarks [Q, and a renormalization group (RG) improved 
gauge action for gluons |T^. The latter choice is motivated by the result that chiral 



symmetry is much better realized with this action than for the plaquette gauge action. We 
may also expect that scaling violation in Bk arising from the gauge action is improved with 
the use of the RG-improved action. 

We examine effects due to finite fifth dimensional size and finite spatial size iVo- in 
detail. Scaling behavior of Bk is studied by adopting (3 = 2.6 and 2.9 corresponding to 
the lattice spacing ~ 2GeV and 3GeV. Matching to the continuum operator is made 
perturbatively at one loop order. Making a constant fit in a for the continuum extrapolation, 
we obtain Bk{^i = 2GeV) = 0.5746(61) (191) as an estimate of the continuum value in the 
MS scheme with naive dimensional regularization (NDR). Here the first error is statistical 
and the second error is an estimated systematic error due to scaling violation and 0(a|jg) 
terms in the renormalization factors. This value is consistent with the Kogut-Susskind result 
quoted above, albeit lying at the lower edge of the one standard deviation error band of the 
latter result. We also report on light quark masses obtained from meson mass measurements 
in our simulation. 

This paper is organized as follows. In section |I| we define the fermion and gluon ac- 
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tions, where we recapitulate the argument for choosing the RG-improved action for gluons. 
Numerical simulations and run parameters are described in section |T|. In section IV we 
discuss the operator matching between the lattice and continuum. Hadron mass results, in 
particular the chiral behavior of pseudo scalar meson mass, are discussed in section 0. Our 
main results for the kaon B parameter are given in section |V1[ Section |V11| is devoted to the 
derivation of light quark mass. We close the paper with a brief summary and comments in 
section IVlli. 



II. ACTION 

We employ Shamir's domain- wall fermion action . Flipping the sign of the Wilson 
term and the domain wall height M, we write 

Sf = - ^i^^s)Ddwfix,s;y,s')^p{y,s') ^ (2-1) 

x,s,y,s' 2; 

Dan^fix, s; y, s') = D\x, y)6s,s' + D\s, s')6,^y + (M - 5)6,J,^s' , (2.2) 
D\x, y) = ^ i [(1 - ^,)U,J,+^,y + (1 + ^,)Ulj,.^ , (2.3) 

r Pl62,s' {s = 1) 

D^{s,s') = I Pl6s+i,s' +PrSs-i,s' (K s < iVs) , (2.4) 

[ Pr6n,-1,s' {s = N5) 

where x, y are four- dimensional space-time coordinates, and s, s' are fifth-dimensional or 
"flavor" indices, bounded as 1 < s, s' < with the free boundary condition at both ends 
(we assume to be even); Pr/l is the projection matrix Pr/l = (1 ± 75)/2, and m/ is 
the bare quark mass. The four-quark operator for our calculation is constructed with the 
4-dimensional quark field defined on the edges of the fifth dimensional space, 

q{x) = PL^ix.l) + PRij{x, N^), 

q{x) = ^(x, N,)Pl + ^{x, 1)Pr. (2.5) 
For the gauge part of the action we employ the following form in 4 dimensions: 



plaquette rectangle chair parallelogram 



9' 



(2.6) 



where the first term represents the standard plaquette action, and the remaining terms are 
six-link loops formed by a 1 x 2 rectangle, a bent 1x2 rectangle (chair) and a 3-dimensional 
parallelogram. The coefficients cq, ■ ■ ■ , C3 satisfy the normalization condition 

Co + 8ci + I6C2 + 8C3 = 1. (2.7) 



The RG-improved action of Iwasaki is defined by setting the parameters to cq 



3.648, ci = — 0.331, C2 = C3 = 0. With this choice of parameters the action is expected 
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to exhibit smooth gauge field fiuctuations approximating those in the continuum hmit bet- 
ter than with the unimproved plaquette action. 

A basic piece of information for our study of with DWQCD is in what range of 
the couphng constant (3 = Q/g^ and domain wall height M DWQCD realizes exact chiral 
symmetry in the limit of infinite fifth dimensional size A^s ^ oo. This point has been 



examined in a number of recent studies [ll4| , p!5| . Investigations using the axial vector Ward- 



Takahashi identity show that a non-zero residual quark mass m^q, which represents chiral 
symmetry breaking, remains even in the limit of infinite fifth dimensional size A^s — > oo if 
the lattice spacing is as coarse as ~1 GeV. 

The chiral property is much improved as the coupling constant is decreased. In the 
range corresponding to ~ 2 GeV, the value of residual quark mass becomes an order 
of magnitude smaller than at a"^ ~ 1 GeV at similar fifth dimensional sizes A^s. For the 
standard plaquette gauge action, it is still not clear whether m^g vanishes exponentially with 
a small decay rate |T3] or remains finite, albeit very small, as A^s — > oo \T^\. In contrast, for 



the RG-improved gauge action, the residual quark mass shows an dependence consistent 
with an exponential decay in up to A^s = 24. Furthermore the magnitude of is an 
order of magnitude smaller than that for the plaquette gauge action. 

We can conclude that chiral symmetry is much better realized with the RG-improved 
gauge action than with the plaquette gauge action. We therefore employ the RG-improved 
gauge action for our investigation of the parameter. 



III. RUN PARAMETERS AND MEASUREMENTS 

Parameters of our simulations and the number of configurations employed are summa- 
rized in Table |. We carry out runs at two values of coupling, P = 2.6 and 2.9, corresponding 
to a lattice spacing = 1.81(4) GeV and 2.81(6) GeV determined from the the p meson 
mass rup = 770MeV. The first value is chosen since chiral symmetry is sufficiently well 



realized \T^\, and the second value is selected to check scaling violation effects. 

For our main runs we use the lattice size A^ x Nf x = 24^ x40xl6at/? = 2.6, and 
32^ x60xl6at/5 = 2.9. These lattices have a reasonably large spatial size of aN„ ~ 2.6 fm 
or 2.3 fm respectively. The choice of = 16 at f3 = 2.6 is based on our previous result 
that the anomalous quark mass is already quite small, m^q = 0.274(42) MeV, for this 
parameter set with the domain wall height M = 1.8. In this paper the domain wall height 
is also taken to be M = 1.8. 

We examine the dependence on the fifth dimensional length at jS = 2.6 for the spatial 
size Ao- = 24 using A^s = 16 and A^s = 32. Since we expect the decay rate in A'5 to become 
larger toward weaker coupling, we only employ A'5 = 16 at /? = 2.9. 

The spatial size dependence is examined at /? = 2.6 varying the spatial size from A^ = 24 
to either A^^. = 16 or 32, which correspond to the physical size of aN„ ~ 1.7 and 3.4 fm. The 
size dependence is also checked at /? = 2.9 by adopting A^. = 24 and 32 (aAg. = 1.7, 2.3 fm). 

We take degenerate quarks in our calculations. The common value of bare quark mass 
is chosen to be m/a = 0.01, 0.02, 0.03, 0.04 at both /3 = 2.6 and 2.9, which covers the range 
mps/mv ~ 0.4 — 0.8. 

Quenched gauge configurations are generated on four-dimensional lattices. A sweep of 
gauge update contains one pseudo-heatbath and four overrelaxation steps. After a ther- 
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malization of 2000 sweeps hadron propagators and 3-point functions necessary to evaluate 
Bk are calculated at every 200th sweep. The gauge configuration on each fifth dimensional 
coordinate s is identical and is fixed to the Coulomb gauge. 

In the course of our simulation we measure the kaon B parameter, 

r. {K\s^,il~^,)ds^,{l-^,)d\K) 

I {K 1^7^75^! 0) (0 \s^,^,d\ K) ^ ■ ' 

and the matrix element divided by the pseudo scalar density, 

R _ {K\s^,{l-^,)dsi,il-l5)d\K) 
^ {K\s^,d\0){0\s^,d\K) 

which should vanish at m^r 0. The s and d quark fields defining these quantities are the 
boundary fields given by (pr5|), and the four-quark and bilinear operators are taken to be 
local in the 4-dimensional space-time. 

The domain-wall quark propagator needed to extract the B parameters above is calcu- 
lated with the conjugate gradient algorithm with an even-odd pre-conditioning. Two quark 
propagators are evaluated for each configuration corresponding to the wall source placed at 
either t = 1 or 40 at /? = 2.6 (t = 4 or 57 at /? = 2.9) in the time direction with the Dirichlet 
boundary condition, while the periodic boundary condition is imposed in the spatial direc- 
tions. The two quark propagators are combined to form the kaon Green function with an 
insertion of the four-quark operator at time slices 1 < t < A^j in a standard manner (see, 
e.g., Ref. §). 

We employ the quark propagators above to also evaluate pseudo scalar and vector meson 
propagators, and extract their masses. These masses are calculated for degenerate quark- 
antiquark pair. The physical point for light quark masses rriud and is calculated by linearly 
fitting the meson masses and my as a function of m/, and using the experimental values 
of mT,/mp and rriKlmp or m^jmp as input. 



IV. OPERATOR MATCHING 

We carry out matching of the lattice and continuum operators at a scale g* = 1/a using 
one-loop perturbation theory [0] and the MS scheme with NDR in the continuum. The 
continuum value at a physical scale e.g., fi = 2 GeV, is obtained via a renormalization group 
running from q* = 1/a to fi. 



5;,(NDR,/i) 



1 - 



47r 



-1 r 



1 - 



4:71 



2/3, 



X 



«Ms(g*) 



-7o/2/3o 



5k (NDR, q* 



(4.1) 



where /3o = 11, /3i = 102, 70 = 4 and 71 = —7 [Tj are the Nj = quenched values for the 
renormalization group coefficients. 

In the domain wall formalism the renormalization factor of an n-quark operator 0„ has 
a generic form 
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Or^(/i) = ZOf «-(l/a), 



lattice / 



(4.2) 

^0) ^o„, (4.3) 

where Wq = 1 — M, and represents the quantum correction to the normahzation factor 
1 — Wq of physical quark fields q,q, and Zo„ is the vertex correction to 0„. In the present 
paper we need the factors Z2, Z^, Za, Zp and Zq^ for the quark wave function, quark 
mass, axial vector current, pseudo scalar density and the four-quark AS* = 2 weak operator. 



Perturbative calculation of these renormalization factors at one loop order is given in Ref. |16 



for the DWQCD system with the standard plaquette gauge action. Here we summarize 
results for the RG-improved gauge action. 

The generic form of the one-loop renormalization factors is given by 

2wn g^C, 



Z^{lia) = 

Z2(yua) = 
Zmil^a) -- 

ZA^liO) = 

Zp{fia) = 
Zoiifia) 



1 + 



JQ 

1 — 

,2 



§ 167r2 



1 + 



1 + 



1 + 



167r2 

g'Cp 

167r2 

g'Cp 

167r2 

g'Cp 

167r2 
-,2 



1 + 



g 



167r2 



logifiaf + Z2{M) 

za{M), 
31og(yua)2 + zp{M) 



(4.4) 

(4.5) 
(4.6) 
(4.7) 
(4.8) 
(4.9) 



where Cp is the second Casimir invariant Cp = 4/3 and the finite part zo„ is a function of 
the domain-wall height M. The difference between the plaquette and the RG action resides 
in the finite part. 

In the first row of Table |I| we list the finite parts of the renormalization factors at 
M = 1.8. The one- loop correction in Z^ is very large for our choice of M because of the 
tadpole factor in Zy^ and division with 1 
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Hence we apply a tadpole improvement 
by explicitly moving the one-loop correction to the domain wall height M from Z^ to Wq 
additively, and by factoring out a tadpole factor u"/^ = p"/8 ^nii^ p tj^g plaquette from 
Zo^- This leads to the rewriting, 

-n/2 



z 



MF 



1 



mfV 
) 



MP 

w 



-n/2 



U 



./2^MF^ 



where 



/7MF 

nrMF 
^2 

^MF 
^•MF 

On 



Wo + 4(1 



u] 



1 - 



g'^CpUi 



Lo=^MF + -g CpUl, 
' ^g^CpUi, 



Zo 



n 



n\W0=wf^ + -^g CpUl. 



(4.10) 

(4.11) 
(4.12) 

(4.13) 
(4.14) 
(4.15) 
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Here ui is the one-loop correction to the tadpole factor u = 1 — g'^CFUi/2 + ■ ■ ■ which has 
the values 

_ J 0.125000 (plaquette action) (am^\ 
~ I 0.052567 (RG improved action) ' ^ ^ ^ 

For the tadpole factor u = P^^^ we use the following value of the plaquette for the RG 
action 

r 0.670632(10) at/3 = 2.6 , . 

~ 1 0.707662(5) at /3 = 2.9 ' ^ ' 



obtained from our main simulations. The domain-wall height is shifted according to ( [4.11| ) 



as 

M-^R^M^^ -{^-^^^^ for/? = 2.6 

In the second and third rows of Table || we list the finite parts of the renormalization factors 
after tadpole improvement. 

A mean-field estimate appropriate for the RG-improved action is used for calculating the 
coupling constant 5'|jg(/i), which is given with the following formula for the quenched case 

1 3 22 

^— = (3.648P - 2.648i?) ^ + -—r \og{fxa) - 0.1006, (4.19) 

where i? is a 1 x 2 rectangular Wilson loop whose value is given as 

r 0.45283(2) at/3 = 2.6 

10.50654(1) at/? = 2.9- ^ ^ 

The gauge coupling at /i = 1/a turns out to be 

, ._r2.2731 at/3 = 2.6 , . 

^Msl VaJ - [ 2.0046 at /5 = 2.9 ' ^ ^ 

For Bk the factor (1 — Wq)'^Z^ cancels out, and the one-loop value is given by the ratio 

7 ( \ l + (-21og(/ia)^ + ^oJ(7V(167r^) , . 9^ ( 9w^,^2 , , \ 
= ^ = (1 + (C,z,),V(16vr^))^ = ' + 16^ (-21og(/.a) + Zb,,) • 

(4.22) 

In Table JTT| we give the finite parts of with and without mean field approximation at 
M = 1.8 together with those for Zsp, 

ZbAH = 1 + (-101og(/ia)2 + ZB,) . (4.23) 

The finite parts zq^ and 2CfZa are very similar in magnitude, albeit individually not very 
small. As a result the finite part zbj^ = zo^—2CfZa for is small, and the renormalization 
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factor for Bk with the tadpole improvement turned out to be very near unity e.g. at the 
matching scale q* = 1/a, 

Z^(a*-l/a)-l^-^^^ at/? = 2.6 

The Z factor at the scale /i = 2GeV obtained with a 2-loop running with Eq. (|4.1|) fL% 
becomes 

7MS. ^p.n _/ 0.979 at/5 = 2.6 
Z^^(/i-2GeV)-|^QQg at/3 = 2.9- 

Meanwhile Zbp is evaluated by setting fi = 2 GeV in Eq. (|4.23| ) 

7MS/ , _ nr^\r\ - / 1-00'^ at /3 = 2.6 
Z5^(/i-2GeV)-|^_^29 at /5 = 2.9 " 

For quark mass the renormalization factor at the matching scale q* = 1/a takes the 
values 

1.173371 at (3 = 2.6 



(4.25) 



(4.26) 



z^i^* = 1/a) = (1 - (^0^^)') (zr) «-^^r (g* = = { 1 



094189 at /3 = 2.9 



. (4.27) 



With a renormalization group running from the scale q* to fi = 2 GeV using the four-loop 
anomalous dimension and beta function |TH1, we have 



(4.28) 



C(«MS(? )/^) 

and the renormalization factor becomes 



S(/i = 2GeV) 



1.155769 at/3 = 2.6 
1.147224 at P = 2.9 ' 



The four-loop running factor c{ajj^{fi) / vr) is given by [|n 



(4.29) 



c{x) = (x)'^"«|l + (7i-/5i7o)a; 
1 



+ 2 



(71 - /9i7o)^ + 72 + /5i%o - /?i7i - /^27o 

^(71 - /3i7o)^ + ^(71 - /3i7o)(72 + A^To - A71 - /?27o) 
6 2 



+ 



where 



l/_ -3_ 
+ (73 - /^i 7o + 2/3i/527o - /337o + Pili - hli - 



28^7 1497^3 
/3o = ll, A = 102, /52 = ^, /33 = ^^ + 3564C(3), 
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7o"^ = 4, 7r 

73™ 



202 



7^ = 1249, 



^^°^°^S1^C(3)- 8800C(5) 



162 



27 



(4.30) 

(4.31) 
(4.32) 
(4.33) 
(4.34) 



with ( the Riemann zeta-function. 

Let us add a comment on the systematic error due to operator matching. Since we have 
used the one-loop renormahzation factor for operator matching, the systematic error should 
include contributions from higher loop corrections. We estimate the magnitude of these 
corrections by changing the matching scale from q* = 1/a to q* = n/a and also adopting a 
different definition for gauge coupling using the plaquette value only given by 

1 3 22 

The gauge coupling at /i = 1/a becomes 

2 n , \ f 1-8839 at P = 2.6 



V. PSEUDO SCALAR AND VECTOR MESON MASSES 

A. Extraction of meson masses 

We extract pseudo scalar and vector meson masses mps and my at each rrif, No- and 
by a single exponential fit of meson propagators. Representative plots of effective mass are 
shown in Figs. |l]and ^. The fitting range chosen from inspection of such plots is 12 < t < 27 
and 6 < t < 16 for pseudo scalar and vector meson mass for all simulations at /5 = 2.6, and 



18 < t < 41 and 10 < t < 26 at /3 = 2.9. In Tables |V[]I3 we list the numerical values of 
rupsa, mya and the ratio at four quark masses irtfa = 0.01,0.02,0.03,0.04 for each set of 
run parameter. The errors given are calculated by a single elimination jackknife procedure. 



B. Chiral extrapolation 



For chiral extrapolation we fit the light hadron masses wip5 and my linearly as a function 
of mja as illustrated in Figs. || and ^ Since pseudo scalar meson mass thus extrapolated 
does not vanish at = 0, we employ a fit of the form 

mpgO^ = Aps{mfa + rriresO.), (5.1) 
mya = Ay + By m fa (5.2) 

and determine the parameters Aps,mres(i for the pseudo scalar meson, and Ay, By for the 
vector meson. The physical point for the bare quark mass parameter mf corresponding to 
physical u and d quark [mj^)-, which are assumed degenerate, and s quark (m^) are fixed 
by the equations 

^ Aps{mfa + mresO) _ m^ 
Ay + BymJ'^a mp 

jAps{m}a/2 + 

Ay + Bym'f'^a mp 

Ay + By m J {(f)) a m^ 

Ay + Bym'^'^a mp 



0.135 
0.77 ' 

0.495 
0.77 ' 

1.0194 
0.77 ' 



(5.3) 

(5.4) 
(5.5) 



10 



where for s quark we employ the kaon (m^) or phi (mj(0)) meson mass as input. We then 
fix the lattice spacing a by setting the vector meson mass at the physical quark mass point 
mj'^ to the experimental value rrip = 770 MeV. Numerical values of lattice spacing and other 
parameters are listed in Table 0. 

In Fig. 1^ we plot results for rripa at the physical point. The values for given (3 are 
reasonably consistent with each other; the variation of results depending on spatial volume 
is mild, and the difference between the fifth dimensional size = 16 and 32 at /5 = 2.6 
on 24^ X 40 lattice is a one-standard deviation effect. In the following analyses we use the 
lattice spacing corresponding to each spatial size and fifth dimensional length. 



C. Chiral property of pseudo scalar meson mass 

We have already mentioned that the pseudo scalar meson mass, if linearly extrapolated, 
does not vanish at = 0. We have also examined alternative fits including either a 
quadratic term, (mja)^, or a quenched chiral logarithm term, mjalog(mja), in addition 
to the hnear term. We have found that these yield almost identical values of the pseudo 
scalar meson mass at = 0. We observe from the results at /? = 2.6 shown in the 
left panel of Fig. |^ that the non-zero pseudo scalar meson mass cannot be explained as 
an effect of finite fifth dimensional lengths, since the data at = 16 (open circles) and 

= 32 (open squares) are consistent within the error down to the smallest quark mass 
rrifa = 0.01. This conclusion is also supported by an analysis of the anomalous quark mass 
m^q defined by the axial Ward-Takahashi identity 0]. This quantity provides a measure of 
chiral symmetry breaking due to a finite N^. It was found that m^g has only a very small 
value of m^g = 0.274(42) MeV for = 16 at l3 = 2.6. For comparison, the magnitude of 
rrires obtained from the linear fit is 2 — 4 MeV as one can see from Table ^. 

Examining the spatial size dependence of results at /3 = 2.6 (left panel of Fig. ^ for 
Ncr = 16, 24 and 32, we observe that the three points are mutually consistent within the 
errors for the heavier quark mass oimja = 0.04, 0.03 and 0.02, but that they show a decrease 
toward larger spatial volumes at our lightest quark mass m^a = 0.01. This indicates that 
the non-zero pseudo scalar meson mass at = in the linear extrapolation refiects a finite 
spatial volume effects in our pseudo scalar meson mass data. 

To make this point explicit, we plot the values of at nif = as a function l/N^a 
in Fig. p. At j3 = 2.6 the results (filled circles) exhibit a decrease as 1/N„a —>■ 0. For 
comparison we plot by open squares results for the Kogut-Susskind quark action, which 
retains U{1) chiral symmetry, obtained at a similar lattice spacing of ~ 2 GeV and 
spatial lattice sizes of N^r ~ 16 — 24 |TP|. A similar magnitude of mj^g in the chiral limit 
between the two quark actions both having chiral symmetry corroborates finite-size effects 
as the origin of non-zero values ^p^. 

The two points for /3 = 2.9 do not show a clear volume dependence. This refiects an 
absence of spatial size dependence at mja = 0.01 — 0.04 observed in the right panel of Fig. |^. 
Quark masses in this range are heavier than those at f3 = 2.6 due to a smaller lattice spacing, 
and hence calculations at smaller values of m/a are needed to expose finite spatial volume 
effects at p = 2.9. 
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VI. B PARAMETERS 



A. Extraction of B parameters 

In Fig. and ^ we show typical data for the ratio of kaon Green functions for Bx and 
Bp defined in and ( |3.2| ) as a function of the temporal site t of the weak operator. The 
values of these quantities at each mj, N^, and are extracted by fitting the plateau with 
a constant. The fitting range, determined by the inspection of plots for the ratio and those 
for the effective pseudo scalar meson mass, is 12 < t < 27 for all simulations at /3 = 2.6 and 
18 < t < 41 at /3 = 2.9. In Tables [rV| - |IX| we list the numerical values of B^ and Bp at four 
quark masses mja = 0.01, 0.02, 0.03, 0.04 for each set of run parameter. 



B. Chiral property for Bp 

We have seen in Sec. |V| that corrections in mp^ due to a finite fifth dimensional size 
is sufficiently small for = 16 for the range of quark mass ruf explored, and that the 
non-zero pseudo scalar meson mass at = is caused by finite spatial size effects. As 
a further check we investigate the chiral property of the matrix element for the four-quark 
operator through Bp, which is expected to vanish linearly at = 0. In Fig. |^ we plot bare 
values of Bp as a function of m^a at /3 = 2.6 and 2.9. Inspecting the results at f3 = 2.6 on 
the left panel of Fig. ^ we observe an agreement for the fifth dimensional size = 16 (open 
circle) and 32 (open square). This shows that = 16 is also large enough for this matrix 
element. 

On the other hand, there is a trend of increase for larger spatial volumes when the quark 
mass goes below rrif = 0.02. Finite spatial size effects appear also in this quantity. Making a 
linear chiral extrapolation, we find a small but negative residual at = 0. Contrary to the 
case of mp5, two alternative fits including a quadratic term, (m^a)^, or a chiral logarithm 
term, rrifa log(mja), give smaller sizes of the intercept at mj = compared to that from 
the linear fit. We find, however, that sizes of the intercept decrease as N^r increases for all 
fits. Therefore we conclude that the non-zero values of Bp at = are a finite-spatial 
size effect, and is not a signal of violation of chiral symmetry. 

The negative sign of the intercept may be understood as follows. Chiral symmetry implies 

(0K^_2m^ 
{0\P\P} mps 

where the bare quantities and P are local axial vector current and pseudo scalar density, 
and Za is the renormalization factor for A^, with which we obtain 

8^ \{0\A^\P)\^ ^ 32m} 32m} 



3 " KO|P|P)P ''3Zlmls ''3ZlAps{mf + 



mr 



This relation is well-satisfied at /3 = 2.6 where Za is non-perturbatively known |^ and is 
reasonably good with the perturbative Za at P = 2.9. Since mres ^ mj in the range of mf 
in our simulation, we approximately obtain 
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32 

Bp ~ Bk^tt^. — {nif - nires), (6.3) 

showing that a positive rrires imphes a negative intercept of Bp. This formula also suggests 
that the large part of the size effect for Bp is caused by that for mj^g. 



C. Bk 

The bare value of Bk is interpolated as a function of rrifa using a formula suggested by 
chiral perturbation theory pT] , 

Bk = B {1 — 3cmja log(mja) + brrifa) . (6.4) 

This interpolation is illustrated in Fig. |10]. The physical value of Bk is obtained at the 
point rrif = mj/2 (solid circles in Fig. ^) which is estimated from the experimental value of 
TTiK/mp. The renormalized values of i?A'(NDR;/i = 2GeV) and related physical quantities 
are collected in Table |X]. 



We plot the renormalized value of Bk as a function of the spatial size in Fig. |rT[ Filled 
circles and triangles are results at /? = 2.6 and 2.9 keeping the same fifth dimensional size 

= 16. At /9 = 2.6 we observe a slight increase of Bk from the spatial size N^a ^ 1.7 fm to 
2.6 fm, but the values beyond the size N^a ^ 2.6 fm are well consistent within the statistical 
error of 1%. This result agrees with that of a previous finite spatial size study with the 
Kogut-Susskind quark action which found finite size effects to be smaller than 0.5% for 
the spatial size N„a>2.2 fm. We conclude that the size of about 2.6 fm {N„ = 24, j3 = 2.6) 
and 2.3 fm {N„ = 32, (3 = 2.9) used in our main runs is sufficient to avoid spatial size effects 
for Bk at a 1% level. 

In Fig. O we plot Bk as a function of the fifth dimensional length iVs on a 24^ x 40 



four-dimensional lattice at P = 2.6. The results at = 32 and = 16 are in agreement 
within the statistical error of 1%. Hence the fifth dimensional size of A^s = 16 is sufficient 
for the calculation of Bk at this accuracy. 



Our final results from the main runs are shown in Fig. [I3| as a function of lattice spacing 
by filled squares. The open symbols and the associated lines represent results from a previous 
calculation with the Kogut-Susskind (staggered) quark action 0, where gauge invariant and 
non-invariant four-quark operators are used. Our result obtained with the domain wall quark 
action and an RG-improved gluon action show a much better scaling behavior; the central 
values of the two points differ by only 1.6% while the Kogut-Susskind results show a 10% 
decrease over the similar range of lattice spacing ^2 — 3 GeV. In order to estimate 
the continuum value, we then make a constant extrapolation BK{a) = Bk-, which yields 
Bk[ii = 2GeV) = 0.5746(61). 

Possible sources of systematic errors in this result are scaling violation ignored in the 
constant fit and higher loop corrections in the renormalization factors. Making an extrapo- 
lation of our data of the form Bk{(i) = Bk + c-o^, based on O(a^) scaling violation expected 



for DWQCD |22,23|, we obtain an estimate of 2.2% for the first error. A simple estimate 
for the second error is provided by the value of aj^(l/a)^ at the finer lattice spacing of 
f3 = 2.9. This yields 2.5% for the second error. This seems to be a reasonable estimate since 
other methods of estimation, either shifting the matching scale from q* = 1/ato q* = vr/a or 
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employing different choices of gauge coupling such as ( [4.35|) , give a small variation of 
Adding the two estimates by quadrature gives a 3.3% systematic error, and we obtain 



5a'(NDR; fx = 2GeV) = 0.5746(61)(191). (6.5) 

as our estimate of the continuum value of in the MS scheme at /i = 2 GeV. 

This value lies at the lower edge of the one-standard deviation error band of the result 
Bk{h = 2GeV) = 0.628(42) obtained with the Kogut-Susskind action 0. We recall that 
the statistical error with the Kogut-Susskind results are at the 0.5 — 1% level. A signifi- 
cantly larger error of 6.7% in the continuum value arises from the continuum extrapolation 
incorporating both the scaling violation and the uncertainty due to the use of one-loop 
renormalization factor. Making a more detailed check of agreement of results from the two 
types of quark actions requires a better control of systematic errors, in particular those due 
to renormalization factors. For this purpose non-perturbative determination of these factors 
for both cases will be necessary. 

The RBC Collaboration carried out a quenched simulation with the domain-wall quark 
action and a plaquette gluon action at /3 = 6.0 (a~^ ~ 2 GeV) on a 16^ X 32 X 16 lattice. 
Employing the method of Ref. to non-perturbatively determine the renormalization fac- 
tors, they reported a value Bxifi = 2GeV) = 0.538(8) [^. This value is 7% smaller than 
our result. A precise comparison, however, would require examination of spatial size and 
scaling violation effects in the RBC result, and of renormalization factors in our result as 
discussed above. 



D. Bk as a function of m^g in the continuum limit 

We have so far discussed the scaling behavior of B^ at the physical quark mass. Our 
data, in fact, allows us to examine the scaling behavior of B^ over a wide range of quark 
mass, and derive the mass dependence of Bk in the continuum limit. 

In order to compare results at different lattice spacings, we employ in physical units 
(GeV^) instead of m^-a. In Fig. |T3| B^ifi = 2GeV) is given as a function of (GeV^) at 



P = 2.6, 2.9 and in the continuum limit. The data are first fitted by 

Bk = B{1 — 3 c m\,g \og{rn?pg) + h rripg) (6.6) 

for each jS and then extrapolated to the continuum by a constant fit. All errors in the 
figure are estimated by a single elimination jackknife procedure, except for fit errors for the 
continuum extrapolation. As seen in the figure, scaling violation is mild up to ^^p^ < 0.8 
GeV^, and the continuum extrapolation is reliable there. This confirms that the small scaling 



violation of the physical Bk observed in Sec. |VI C| is not an accidental one at mps = rriK 



but it holds over a wide range of the pseudo scalar meson mass. 



In Table pCII| , values of Bk in the continuum limit, which are also fitted by the same 
form ( |6.6| ), are given for 0.02 < mp^ < 1.0 (GeV^) with errors. Fitted parameters B, b 
and c are also given in the table, together with the reproduced values. From this result 
in the continuum limit one can see that the contribution from higher order terms in chiral 
perturbation theory [b and c) is non-negligible and becomes as large as 40% of the leading 
order contribution {B) at rups = rriK- 
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VII. LIGHT QUARK MASSES 



We attempt a determination of light quark masses niud = ijnu + rndj/^ and using 
our meson mass data. There is a difficulty associated with a non-zero pseudo scalar meson 
mass at = due to finite spatial sizes, which is represented by m^es in the linear chiral 
formula ( |5.1| ). This causes systematic uncertainties in the results for quark masses, which is 
quite sizable for light u and d quarks. 

In order to examine this problem, we calculate the physical quark masses in two ways 
which differ in the choice of origin for bare quark mass. In the first method we take = 
as the origin, and write 



m^d = Z,m)\ (7.1) 



ms = ZM-mf). (7.2) 



Here m'^'^ and rrij are the bare quark mass rUf for the physical point of pion and kaon 
determined by (^.3|) and (5^). The subtraction of mj'^ in the second equation is to take into 



account the contribution of u-d quark in the kaon mass, mx oc rriud + ^s, and Zq denotes 
the renormalization factor to match the bare lattice value to that in the continuum in the 
MS scheme with NDR at yU = 2 GeV as discussed in Sec. |V]. 

In the second case we take the point rrif = —rrires, where pseudo scalar meson mass 
vanishes, as the origin. The formula then reads 

mud = Zq{m'f + nires), (7.3) 
ms = Zq{m)-mf + mres)- (7.4) 

On the other hand the strange quark mass with the phi meson mass as input is given directly 



as 



rus = Zq'ms{(j)). (7.5) 



The results of these calculations are listed in Table XIII. 



In Fig. |15| we plot the u-d quark mass calculated in the two ways above as a function of 
spatial size aNcr in physical units. Two features are quite evident from this figure: (i) There 
is little dependence on the fifth dimensional size N^. Hence = 16 is sufficient to avoid 
effects of chiral symmetry breaking at our range of lattice spacings. (ii) Effects of finite 
spatial size, by contrast, are quite significant if rrires is ignored, even yielding a negative 



value for rriud for small spatial sizes (left panel of Fig. |T5]). The values calculated including 
rrires, on the other hand, are much more stable as a function of aN^r (right panel). 

In order to understand the second point, we note that rriresO' depends strongly on the 
volume while the slope Aps is almost volume independent. Using (|5.3| ) and the corresponding 
one at A^o- = oo given by 



^ApsjN^ = oo) ■ mf{N^ = ^ ^ 

Ay + Bvm'j'^a rrip ' 



(7.6) 



and neglecting a small volume dependence of the denominator Ay + Bvm'j'^a, we observe 
that the following formula holds: 
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mf - mf{N, = oo) _ m^es Aps - Aps{N, = oo) 

mf{N^ = oo) m'f{N^ = 00) Aps ' ^ ' ^ 

Since the magnitude m^es ~ 2 — 3 MeV for our spatial size of aN„ 2.5 fm is comparable to 
the actual u-d quark mass, the first term is 0{1) and becomes the main contribution to the 
size effect, while the second term, representing finite size effect in the slope Aps, is found to 
be much smaller. Hence including rrires removes a dominant part of finite size effects in u-d 
quark mass. In view of this situation we take the values including the best estimate 

from our present data for rriud- 

The scaling behavior of rriud is plotted in Fig. ^ by filled squares. Making a constant fit 
to the two values, we find 

m^(2MeV) = 3.764(81) (215)MeV, (7.8) 

where the first error is statistical and the second due to scaling violation and 0{a'^) system- 
atic errors estimated in the same way as for Bk (see Sec. [VI Q ) and added in quadrature. 

For the heavier s quark, effects of less significant, and those of are within the 

statistical error, as one can see in Table pCIIl]. In parallel with u-d quark mass we take the 



results including m^es as our best estimate. The two values from our main runs at /? = 2.6 
and 2.9 are plotted by filled squares in Fig. {K input) and (0 input). Fitting with a 
constant and making estimation of systematic errors as for m^^ we obtain 

In Figs. 0, |1^ and |18] open symbols show results obtained with the conventional 4- 
dimensional quark actions; circles and squares for the Wilson action with the plaquette gluon 
action p^], diamonds and down triangles for the clover action with the RG-improved gluon 



action as used in the present work ||25|, and right triangles for the Kogut-Susskind quark 



action with the plaquette gluon action PB[. The first two cases use one-loop renormalization 



factors, while the Kogut-Susskind results are based on a non-perturbative value calculated 
in the RI scheme. The values estimated in the continuum limit in each of these studies are 



plotted at a = and are summarized in Table pCIV . 

Compared to the values obtained with the 4-dimensional quark actions, our results with 
the domain-wall action are somewhat small both for u-d and s quark. As with the case of 
Bk, a more precise examination of the issue of agreement of the continuum value requires a 
non-perturbative determination of the renormalization factors for our combination of quark 
and gluon actions. 



We note that a recent result = 110 (2) (22) MeV [^] with K input using domain wall 
fermions and non-perturbative renormalization factor but with the plaquette gauge action 
at /3 = 6.0 is consistent with ours, within the 20% systematic error quoted which includes 
that associated with the conversion from the RI scheme to the MS scheme. 



VIII. CONCLUSIONS 

In this article we have presented our investigation of quenched calculation of the kaon B 
parameter B^ with domain-wall QCD. 
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In order to make full use of the good chiral property of this system, we employed a 
renormalization-group improved gluon action and carried out simulations at a~^>2 GeV. 
According to our previous study [Q, the magnitude of chiral symmetry breaking due to 
finite fifth dimensional size N^,, if measured in terms of residual quark mass m^q, is less than 
1 MeV for A'"5>10 at such lattice spacings. An explicit examination of the dependence 
of Bk has shown that such effect is less than 1% for > 16 at a~^>2 GeV. 

We have also found that spatial size effects are less than 1% for the physical spatial sizes 
aN„>2.5 fm, confirming the finding of a previous study with the Kogut-Susskind quark 
action. Furthermore, scaling violation turned out to be very small, being less than 2% 
between a^^ ^ 2 GeV and 3 GeV. 

These results show that DWQCD, albeit computer time consuming by a factor O^N^) 
compared to conventional lattice QCD simulations, provides a very good framework for a 
precision determination of B^- An important ingredient toward this goal, which was not 
available for the present study, is the value of the renormalization factors precise to the level 
of one percent. Results using the RI scheme have been reported for the plaquette action by 
the RBC Collaboration , and an attempt employing the Schrodinger functional technique 
is in progress pO[ . Hopefully progress in these calculations will allow us to report results for 
Bk in the continuum limit with a total error of at most a few percent in quenched QCD in 
the near future. 

We have also examined the possibility of calculating light quark masses in DWQCD. We 
find good scaling behavior, and the values estimated for the continuum limit are in reasonable 
agreement, albeit somewhat small, with those of 4-dimensional simulations. Further progress 
toward precision determination of light quark masses also requires that of renormalization 
factors at the few percent level. 
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TABLES 





2.6 


2.6 


2.6 


2.6 


2.9 


2.9 


a-^GeV) 


1.875(56) 


1.807(37) 


1.758(51) 


1.847(43) 


2.869(68) 


2.807(55) 


Nt 


40 


40 


40 


40 


60 


60 


Na 


16 


24 


24 


32 


24 


32 


N5 


16 


16 


32 


16 


16 


16 


N^a (fm) 


1.7 


2.6 


2.6 


3.4 


1.7 


2.3 


# conf. 


122 


76 


50 


25 


76 


50 



TABLE I. Simulation parameters together with the number of configurations analyzed shown 
in bold letters. 



M 


Z2 






ZA 


Zp 


ZO4, 


1.8 


-3.824 


13.148 


-25.1295 


-9.190 


-13.148 


-23.868 


M 




~MF 


~MF 


~MF 


~MF 
Zp 


~MF 


1.4198 


0.651 


6.044 


-7.92355 


-4.692 


-6.044 


-13.612 


1.4687 


0.632 


6.319 


-7.95874 


-4.714 


-6.319 


-13.500 



TABLE II. Finite parts of the renormalization factors with RG improved gauge action. The 
mean field approximation is used for the factors at M = 1.4198 and 1.4687. Errors from the 
numerical integration are in the last written digit. 



M 


ZBk 


ZBp 


~MF 


1.8 


0.64 


11.19 




1.4198 


-1.10 




2.51 


1.4687 


-0.93 




3.35 



TABLE III. Finite parts of the renormalization factors zbj^ and ZBp with RG improved gauge 
action. The mean field approximation is used for ZBp at M = 1.4198 and 1.4687, while the effect 
of the approximation on zbj^ is just to shift the domain- wall height. Errors from the numerical 
integration are in the last written digit. 
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(3 = 2.6 on 16^ x 40 x 16 lattice 


rrif 


Bk 


Bp 


mps 


my 


mps/mv 


0.01 


0.487(19) 


0.0225(10) 


0.1902(29) 


0.438(12) 


0.434(13) 


0.02 


0.566(14) 


0.0569(14) 


0.2554(23) 


0.4639(72) 


0.5505(97) 


0.03 


0.615(11) 


0.0945(17) 


0.3083(21) 


0.4904(54) 


0.6287(80) 


0.04 


0.6486(85) 


0.1327(19) 


0.3548(20) 


0.5168(43) 


0.6865(69) 



TABLE IV. Data for Bk, Bp, mps, my and mps/mv at each quark mass nif at /? = 2.6 on 
16^ X 40 X 16 lattice. 





P = 2.6 on 243 X 40 X 16 lattice 


nif 


Bk 


Bp 


mps 


my 


mps/mv 


0.01 


0.5196(69) 


0.02487(38) 


0.1883(13) 


0.4522(83) 


0.4164(80) 


0.02 


0.5769(47) 


0.05777(55) 


0.2567(12) 


0.4683(51) 


0.5481(65) 


0.03 


0.6155(42) 


0.09367(70) 


0.3102(11) 


0.4900(36) 


0.6331(52) 


0.04 


0.6453(39) 


0.13092(83) 


0.35647(97) 


0.5141(29) 


0.6934(42) 



TABLE V. Data for Bk, Bp, mps, my and mps/mv at each quark mass nif at P = 2.6 on 
24^ X 40 X 16 lattice. 





P = 2.6 on 243 X 40 X 32 lattice 


nif 


Bk 


Bp 


mps 


my 


mps/mv 


0.01 


0.5225(89) 


0.02444(47) 


0.1874(16) 


0.463(12) 


0.405(11) 


0.02 


0.5780(63) 


0.05751(71) 


0.2558(14) 


0.4760(73) 


0.5374(89) 


0.03 


0.6152(54) 


0.09329(90) 


0.3093(12) 


0.4956(50) 


0.6242(67) 


0.04 


0.6437(49) 


0.1302(10) 


0.3556(11) 


0.5182(37) 


0.6862(53) 



TABLE VI. Data for Bk, Bp, mps, my and mps/mv at each quark mass nif at /? = 2.6 on 
24'^ X 40 X 32 lattice. 



21 







(5 = 2.6 on 323 


X 40 X 16 lattice 






rrif 


Bk 


Bp 


mps 


my 


mps/mv 


0.01 


0.5260(47) 


0.02594(33) 


0.1841(13) 


0.4424(92) 


0.4161(86) 


0.02 


0.5766(40) 


0.05873(51) 


0.2535(12) 


0.4608(59) 


0.5501(77) 


0.03 


0.6124(39) 


0.09460(72) 


0.3076(12) 


0.4847(44) 


0.6347(66) 


0.04 


0.6424(35) 


0.13196(88) 


0.3544(11) 


0.5097(36) 


0.6954(57) 



TABLE VII. Data for Bk, Bp, mps, my and mps/mv at each quark mass m/ at /3 = 2.6 on 
32^ X 40 X 16 lattice. 





P = 2.9 on 243 X 60 X 16 lattice 


nif 


Bk 


Bp 


mps 


mv 


mps/mv 


0.01 


0.512(18) 


0.0340(14) 


0.1444(20) 


0.2915(59) 


0.496(13) 


0.02 


0.591(11) 


0.0802(17) 


0.2007(17) 


0.3172(40) 


0.6328(96) 


0.03 


0.6403(84) 


0.1295(21) 


0.2460(15) 


0.3425(32) 


0.7182(78) 


0.04 


0.6758(75) 


0.1793(24) 


0.2856(14) 


0.3681(27) 


0.7760(65) 



TABLE VIII. Data for Bk, Bp, mps, my and mpslray at each quark mass mj at /? = 2.9 
on 24^ X 60 X 16 lattice. 





/3 = 2.9 on 32^ x 60 x 16 lattice 


my 


Bk 


Bp 


mps 


my 


mps/my 


0.01 


0.5318(72) 


0.03412(64) 


0.1459(12) 


0.2984(54) 


0.489(10) 


0.02 


0.5922(58) 


0.07885(97) 


0.2022(12) 


0.3209(31) 


0.6302(78) 


0.03 


0.6345(51) 


0.1271(12) 


0.2470(12) 


0.3451(23) 


0.7157(65) 


0.04 


0.6669(46) 


0.1761(14) 


0.2863(11) 


0.3700(19) 


0.7738(55) 



TABLE IX. Data for Bk, Bp, mps, my and mps/my at each quark mass m/ at /? = 2.9 on 
32^ X 60 X 16 lattice. 
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2.6 


2.6 


2.6 


2.6 


2.9 


2.9 




40 


40 


40 


40 


60 


60 


N 


16 


24 


24 


32 


24 


32 




16 


16 


32 


16 


16 


16 


a"^(GeV) 


1.875(56) 


1.807(37) 


1.758(51) 


1.847(43) 


2.869(68) 


2.807(55) 


mv{mf = 0)a 


0.411(12) 


0.4261(87) 


0.438(13) 


0.4169(97) 


0.2684(64) 


0.2743(54) 


mlsirrif = 0)a^ 


0.0060(12) 


0.00490(53) 


0.00467(69) 


0.00330(52) 


0.00038(62) 


0.00099(34) 




0.00201(42) 


0.00161(18) 


0.00153(23) 


0.00108(17) 


0.00019(31) 


0.00049(17) 


m^es (MeV) 


3.77(78) 


2.90(32) 


2.70(42) 


1.99(32) 


0.54(89) 


1.38(48) 




-0.00027(41) 


0.00022(19) 


0.00040(26) 


0.00067(18) 


0.00091(32) 


0.00066(18) 


m)a/2{K) 


0.0216(15) 


0.0233(10) 


0.0248(15) 


0.0227(11) 


0.01475(82) 


0.01515(70) 


mja{(f)) 


0.0503(62) 


0.0632(64) 


0.071(11) 


0.0583(59) 


0.0350(25) 


0.0373(25) 



TABLE X. Results of meson mass fits. 



/5 


2.6 


2.6 


2.6 


2.6 


2.9 


2.9 


Nt 


40 


40 


40 


40 


60 


60 




16 


24 


24 


32 


24 


32 


N5 


16 


16 


32 


16 


16 


16 


a-^GeV) 


1.875(56) 


1.807(37) 


1.758(51) 


1.847(43) 


2.869(68) 


2.807(55) 


bare B parameters 












Bk 


0.575(14) 


0.5908(57) 


0.5975(77) 


0.5871(60) 


0.554(14) 


0.5655(69) 


renormalized B parameters 


(MS scheme with NDR at = 


2 GeV) 






BKiq* = 1/a) 


0.564(14) 


0.5782(55) 


0.5839(75) 


0.5753(58) 


0.558(14) 


0.5690(70) 


BK{q* = vr/a) 


0.570(14) 


0.5844(56) 


0.5901(76) 


0.5815(59) 


0.563(15) 


0.5741(70) 


BK{g = gp) 


0.566(14) 


0.5803(56) 


0.5862(75) 


0.5773(59) 


0.557(14) 


0.5684(70) 



TABLE XL Results for B parameters together with those of relevant c 
last row denotes the use of an alternative definition of the coupling (4.35) 



5fp in the 
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Ucli dlllC. LCi o 


Bk 


error 


J-'K 


R 






41 1 80 


u 






n 71 1 1 n 


C 






D 9n7'^1 


PS / 


continuum extrapolation 


lecoiioLi ucLion uy Liie iii 


u.uzu 


0.4318 


0.0100 


n A VT7 


n 1 nn 

U.iUU 


0.4994 


0.0058 


U.OUUi 


u.zuu 


0.5544 


0.0044 


u.oozo 


n "inn 


0.5937 


0.0043 


U.oyzz 


0.400 


0.6239 


0.0041 


0.6228 


0.500 


0.6471 


0.0038 


0.6470 


0.600 


0.6643 


0.0038 


0.6660 


0.700 


0.6780 


0.0043 


0.6807 


0.800 


0.6907 


0.0051 


0.6918 


0.900 


0.7016 


0.0060 


0.6996 


1.000 


0.7097 


0.0075 


0.7046 



TABLE XII. Parameters for the fit of Bk in the continuum limit by Eq. ( |6.6D , and Bk as a 
function of rr?pg in the continuum limit, together with the reconstruction from the fit. 





2.6 2.6 2.6 2.6 


2.9 2.9 


Nt 


40 40 40 40 
16 24 24 32 
16 16 32 16 


60 60 
24 32 
16 16 


a-i(GeV) 


1.875(56) 1.807(37) 1.758(51) 1.847(43) 


2.869(68) 2.807(55) 


renormalized quar 
m„,rf(MeV) 
nisiMeV) 


-c masses ignoring mrcs 

-0.60(90) 0.47(39) 0.82(52) 1.44(37) 
94.9(33) 96.8(22) 99.8(31) 95.9(24) 


3.0(10) 2.12(57) 
94.4(27) 95.4(23) 


renormalized quar 


< masses including rrircs 



mu,d + mres(MeV) 
rus + mres(MeV) 


3.79(12) 3.821(79) 
99.3(33) 100.2(21) 


3.92(11) 
102.9(30) 


3.748(92) 
98.3(24) 


3.625(93) 
95.0(24) 


3.701(83) 
97.0(22) 


renormalized s quark mass with (j) input 

ms{MeV) 110.(10) 132.(11) 


144.(19) 


124.8(98) 


115.6(58) 


120.0(57) 


TABLE XIII. 


Results for quark masses. 


Renormalizec 


values are in the MS 


scheme at 



H = 2 GeV. 



24 



ref. quark action gluon action 


mff{2GeV) mf^{2GeV) 

K input input 


this work DW RG-improved 


3.764(81)(215) MeV 98.7(2.1)(5.6) MeV 122.6(6.8)(13) MeV 


[||] Wilson plaquette 

clover RG-improved 
(ID KS plaquette 


4.57(18) MeV 116 (3) MeV 144 (6) MeV 
4.36t|j:}^ MeV llOl^ MeV 132^^ MeV 
4.23(29) MeV 106 (7) MeV 129 (12) MeV 



TABLE XIV. Results for light quark masses as compared with previous studies. One-loop 
approximation to the renormalization factors are employed except for those with the KS fermion 
action in the last row. 
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FIG. 1. Effective pseudo scalar meson mass as a function of temporal distance t &t (3 = 2.6 on 
a 24^ X 40 X 16 lattice (left) and at /? = 2.9 on a 32^ x 60 x 16 lattice (right). Lines show constant 
fit over the fitted range. 
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FIG. 2. Effective vector meson mass as a function of temporal distance t at P = 2.6 on a 
24^ X 40 X 16 lattice (left) and at /? = 2.9 on a 32^ x 60 x 16 lattice (right). Lines show constant 
fit over the fitted range. 
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FIG. 3. Pseudo scalar meson mass squared as a function of bare quark mass at /? = 2.6 
(left) and at /? = 2.9 (right). Lines show linear fits to main runs. The data except for the main 
run are shifted in m/. 
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FIG. 4. Vector meson mass as a function of bare quark mass m/ at /? = 2.6 (left) and at j3 = 2.9 
(right). Lines show linear fits to main runs. The data except for the main run are shifted in ruf. 
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FIG. 5. Rho meson mass as a function of spatial lattice size iV^. at /? = 2.6 (left) and at /? = 2.9 
(right). Filled symbols represent data at fifth dimensional length = 16 and an open circle 
represents that at iVs = 32. 
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FIG. 6. Pseudo scalar meson mass at mj = as a function of spatial lattice size N„. Filled 
symbols represent data at fifth dimensional length = 16 and an open circle represents that at 
A^5 = 32. The results of the KS fermion at similar volume size |19] is also plotted with open squares 
for comparison. 
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FIG. 7. Ratio of weak matrix element with vacuum saturation ( p. ID as a function of temporal 
distance t ai {3 = 2.6 on a 24^ x 40 x 16 lattice (left) and at /3 = 2.9 on a 32^ x 60 x 16 lattice 
(right). Lines show constant fit over the fitted range. 
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FIG. 8. Ratio of weak matrix element with pseudo scalar density ( |3.2| ) as a function of temporal 
distance t at l3 = 2.6 on a 24^ x 40 x 16 lattice (left) and at /? = 2.9 on a 32^ x 60 x 16 lattice 
(right). Lines show constant fit over the fitted range. 




30 



0.8 



0.7 



0.6 



• B,^ at physical quark mass 



0.5 



0.4 



0.8 



0.7 



0.6 



• B,^ at physical quark mass 



CO 



0.5 



0.4 



0.3 



0.02 



0.04 



0.3 



0.02 



0.04 



FIG. 10. Bare Bk interpolated as a function of nifa at P = 2.6 for a 24^ x 40 x 16 lattice (left) 
and at P = 2.9 on a 323 x 60 x 16 lattice (right). 
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FIG. 11. Renormalized Bk as a function of spatial size. 
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FIG. 12. Renormalized Bk as a function of fifth dimensional length N^. 
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FIG. 13. Scaling behavior of renormalized B^ifJ- = 2GeV). Previous results with the KS action 
are also shown with open symbols. 



32 



0.8 



0.7 




p=2.6 

- p=2.9 
o |3=2.6 
o p=2.9 

- a=0 

• physical point 



" 0.2 0.4 0.6 0.8 1 

mps (GeV') 

FIG. 14. Rcnormalizcd Bk{i-I' = 2GeV) as a function of mpg, at /3 = 2.6, 2.9 and cx) (continuum 
limit), where errors are shown only in the continuum hmit. Open symbols represent data obtained 
in our simulations, while the solid circle gives the value of Bk at the physical point in the continuum 
limit, with the statistical error(solid) and the total error (dotted). 
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FIG. 15. Renormalized u-d quark mass as a function of spatial size. Results with (right) and 
without TOres added (left) are shown. 
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FIG. 16. Scaling behavior of u-d quark mass, calculated with nires added, compared with 
those from 4-dimensional quark action: Wilson action (Std) |^4|, clover-improved action (Imp) 
|25|, and Kogut-Susskind (KS) ||2^ actions. VWI and AWI represent vector and axial- vector 
Ward-Takahashi identity masses, respectively, p for the KS fermion represents the matching scale 
of the RI scheme in units of GeV. 



140 



130 



^ 120 

> 

110 

> 

(3 

cvj 100 
II 

E 90 



80 



70 



0.0 



0.5 



1.0 



■ DWF 

o VWI, qStd 
□ AWI, qStd 
OVWI, qImp 
VAWI, qImp 
> KS, p=2.6 



1.5 



a [GeV ] 

FIG. 17. Scaling behavior of the strange quark mass with K input, calculated with mi-cs added, 
compared with those from 4-dimensional quark action: Wilson action (Std) |^^, clover-improved 
action (Imp) |2|], and Kogut-Susskind (KS) |^ actions. VWI and AWI represent vector and 
axial-vector Ward-Takahashi identity masses, respectively, p for the KS fermion represents the 
matching scale of the RI scheme in units of GeV. 
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FIG. 18. Scaling behavior of the strange quark mass with (j) input, compared with 
those from 4-dimensional quark action: Wilson action (Std) [^4|, clover-improved action (Imp) 
1 25 1, and Kogut-Susskind (KS) ||2^ actions. VWI and AWI represent vector and axial- vector 
Ward-Takahashi identity masses, respectively, p for the KS fermion represents the matching scale 
of the RI scheme in units of GeV. 
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